\(\int \frac {(e+f x) \sin ^3(c+d x)}{a+a \sin (c+d x)} \, dx\) [193]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [C] (verified)
   Fricas [A] (verification not implemented)
   Sympy [B] (verification not implemented)
   Maxima [F(-2)]
   Giac [B] (verification not implemented)
   Mupad [B] (verification not implemented)

Optimal result

Integrand size = 26, antiderivative size = 158 \[ \int \frac {(e+f x) \sin ^3(c+d x)}{a+a \sin (c+d x)} \, dx=\frac {3 e x}{2 a}+\frac {3 f x^2}{4 a}+\frac {(e+f x) \cos (c+d x)}{a d}+\frac {(e+f x) \cot \left (\frac {c}{2}+\frac {\pi }{4}+\frac {d x}{2}\right )}{a d}-\frac {2 f \log \left (\sin \left (\frac {c}{2}+\frac {\pi }{4}+\frac {d x}{2}\right )\right )}{a d^2}-\frac {f \sin (c+d x)}{a d^2}-\frac {(e+f x) \cos (c+d x) \sin (c+d x)}{2 a d}+\frac {f \sin ^2(c+d x)}{4 a d^2} \]

[Out]

3/2*e*x/a+3/4*f*x^2/a+(f*x+e)*cos(d*x+c)/a/d+(f*x+e)*cot(1/2*c+1/4*Pi+1/2*d*x)/a/d-2*f*ln(sin(1/2*c+1/4*Pi+1/2
*d*x))/a/d^2-f*sin(d*x+c)/a/d^2-1/2*(f*x+e)*cos(d*x+c)*sin(d*x+c)/a/d+1/4*f*sin(d*x+c)^2/a/d^2

Rubi [A] (verified)

Time = 0.15 (sec) , antiderivative size = 158, normalized size of antiderivative = 1.00, number of steps used = 11, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.269, Rules used = {4611, 3391, 3377, 2717, 3399, 4269, 3556} \[ \int \frac {(e+f x) \sin ^3(c+d x)}{a+a \sin (c+d x)} \, dx=\frac {f \sin ^2(c+d x)}{4 a d^2}-\frac {f \sin (c+d x)}{a d^2}-\frac {2 f \log \left (\sin \left (\frac {c}{2}+\frac {d x}{2}+\frac {\pi }{4}\right )\right )}{a d^2}+\frac {(e+f x) \cos (c+d x)}{a d}+\frac {(e+f x) \cot \left (\frac {c}{2}+\frac {d x}{2}+\frac {\pi }{4}\right )}{a d}-\frac {(e+f x) \sin (c+d x) \cos (c+d x)}{2 a d}+\frac {3 e x}{2 a}+\frac {3 f x^2}{4 a} \]

[In]

Int[((e + f*x)*Sin[c + d*x]^3)/(a + a*Sin[c + d*x]),x]

[Out]

(3*e*x)/(2*a) + (3*f*x^2)/(4*a) + ((e + f*x)*Cos[c + d*x])/(a*d) + ((e + f*x)*Cot[c/2 + Pi/4 + (d*x)/2])/(a*d)
 - (2*f*Log[Sin[c/2 + Pi/4 + (d*x)/2]])/(a*d^2) - (f*Sin[c + d*x])/(a*d^2) - ((e + f*x)*Cos[c + d*x]*Sin[c + d
*x])/(2*a*d) + (f*Sin[c + d*x]^2)/(4*a*d^2)

Rule 2717

Int[sin[Pi/2 + (c_.) + (d_.)*(x_)], x_Symbol] :> Simp[Sin[c + d*x]/d, x] /; FreeQ[{c, d}, x]

Rule 3377

Int[((c_.) + (d_.)*(x_))^(m_.)*sin[(e_.) + (f_.)*(x_)], x_Symbol] :> Simp[(-(c + d*x)^m)*(Cos[e + f*x]/f), x]
+ Dist[d*(m/f), Int[(c + d*x)^(m - 1)*Cos[e + f*x], x], x] /; FreeQ[{c, d, e, f}, x] && GtQ[m, 0]

Rule 3391

Int[((c_.) + (d_.)*(x_))*((b_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Simp[d*((b*Sin[e + f*x])^n/(f^2*n^
2)), x] + (Dist[b^2*((n - 1)/n), Int[(c + d*x)*(b*Sin[e + f*x])^(n - 2), x], x] - Simp[b*(c + d*x)*Cos[e + f*x
]*((b*Sin[e + f*x])^(n - 1)/(f*n)), x]) /; FreeQ[{b, c, d, e, f}, x] && GtQ[n, 1]

Rule 3399

Int[((c_.) + (d_.)*(x_))^(m_.)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(n_.), x_Symbol] :> Dist[(2*a)^n, Int[(c
 + d*x)^m*Sin[(1/2)*(e + Pi*(a/(2*b))) + f*(x/2)]^(2*n), x], x] /; FreeQ[{a, b, c, d, e, f, m}, x] && EqQ[a^2
- b^2, 0] && IntegerQ[n] && (GtQ[n, 0] || IGtQ[m, 0])

Rule 3556

Int[tan[(c_.) + (d_.)*(x_)], x_Symbol] :> Simp[-Log[RemoveContent[Cos[c + d*x], x]]/d, x] /; FreeQ[{c, d}, x]

Rule 4269

Int[csc[(e_.) + (f_.)*(x_)]^2*((c_.) + (d_.)*(x_))^(m_.), x_Symbol] :> Simp[(-(c + d*x)^m)*(Cot[e + f*x]/f), x
] + Dist[d*(m/f), Int[(c + d*x)^(m - 1)*Cot[e + f*x], x], x] /; FreeQ[{c, d, e, f}, x] && GtQ[m, 0]

Rule 4611

Int[(((e_.) + (f_.)*(x_))^(m_.)*Sin[(c_.) + (d_.)*(x_)]^(n_.))/((a_) + (b_.)*Sin[(c_.) + (d_.)*(x_)]), x_Symbo
l] :> Dist[1/b, Int[(e + f*x)^m*Sin[c + d*x]^(n - 1), x], x] - Dist[a/b, Int[(e + f*x)^m*(Sin[c + d*x]^(n - 1)
/(a + b*Sin[c + d*x])), x], x] /; FreeQ[{a, b, c, d, e, f}, x] && IGtQ[m, 0] && IGtQ[n, 0]

Rubi steps \begin{align*} \text {integral}& = \frac {\int (e+f x) \sin ^2(c+d x) \, dx}{a}-\int \frac {(e+f x) \sin ^2(c+d x)}{a+a \sin (c+d x)} \, dx \\ & = -\frac {(e+f x) \cos (c+d x) \sin (c+d x)}{2 a d}+\frac {f \sin ^2(c+d x)}{4 a d^2}+\frac {\int (e+f x) \, dx}{2 a}-\frac {\int (e+f x) \sin (c+d x) \, dx}{a}+\int \frac {(e+f x) \sin (c+d x)}{a+a \sin (c+d x)} \, dx \\ & = \frac {e x}{2 a}+\frac {f x^2}{4 a}+\frac {(e+f x) \cos (c+d x)}{a d}-\frac {(e+f x) \cos (c+d x) \sin (c+d x)}{2 a d}+\frac {f \sin ^2(c+d x)}{4 a d^2}+\frac {\int (e+f x) \, dx}{a}-\frac {f \int \cos (c+d x) \, dx}{a d}-\int \frac {e+f x}{a+a \sin (c+d x)} \, dx \\ & = \frac {3 e x}{2 a}+\frac {3 f x^2}{4 a}+\frac {(e+f x) \cos (c+d x)}{a d}-\frac {f \sin (c+d x)}{a d^2}-\frac {(e+f x) \cos (c+d x) \sin (c+d x)}{2 a d}+\frac {f \sin ^2(c+d x)}{4 a d^2}-\frac {\int (e+f x) \csc ^2\left (\frac {1}{2} \left (c+\frac {\pi }{2}\right )+\frac {d x}{2}\right ) \, dx}{2 a} \\ & = \frac {3 e x}{2 a}+\frac {3 f x^2}{4 a}+\frac {(e+f x) \cos (c+d x)}{a d}+\frac {(e+f x) \cot \left (\frac {c}{2}+\frac {\pi }{4}+\frac {d x}{2}\right )}{a d}-\frac {f \sin (c+d x)}{a d^2}-\frac {(e+f x) \cos (c+d x) \sin (c+d x)}{2 a d}+\frac {f \sin ^2(c+d x)}{4 a d^2}-\frac {f \int \cot \left (\frac {c}{2}+\frac {\pi }{4}+\frac {d x}{2}\right ) \, dx}{a d} \\ & = \frac {3 e x}{2 a}+\frac {3 f x^2}{4 a}+\frac {(e+f x) \cos (c+d x)}{a d}+\frac {(e+f x) \cot \left (\frac {c}{2}+\frac {\pi }{4}+\frac {d x}{2}\right )}{a d}-\frac {2 f \log \left (\sin \left (\frac {c}{2}+\frac {\pi }{4}+\frac {d x}{2}\right )\right )}{a d^2}-\frac {f \sin (c+d x)}{a d^2}-\frac {(e+f x) \cos (c+d x) \sin (c+d x)}{2 a d}+\frac {f \sin ^2(c+d x)}{4 a d^2} \\ \end{align*}

Mathematica [A] (verified)

Time = 6.79 (sec) , antiderivative size = 298, normalized size of antiderivative = 1.89 \[ \int \frac {(e+f x) \sin ^3(c+d x)}{a+a \sin (c+d x)} \, dx=\frac {\left (\cos \left (\frac {1}{2} (c+d x)\right )+\sin \left (\frac {1}{2} (c+d x)\right )\right ) \left (\sin \left (\frac {1}{2} (c+d x)\right ) \left (8 d (e+f x) \cos (c+d x)-f \cos (2 (c+d x))+2 \left (-8 d e+6 c d e+4 c f-3 c^2 f+6 d^2 e x-4 d f x+3 d^2 f x^2-8 f \log \left (\cos \left (\frac {1}{2} (c+d x)\right )+\sin \left (\frac {1}{2} (c+d x)\right )\right )-4 f \sin (c+d x)-d (e+f x) \sin (2 (c+d x))\right )\right )+\cos \left (\frac {1}{2} (c+d x)\right ) \left (8 d (e+f x) \cos (c+d x)-f \cos (2 (c+d x))+2 \left (6 c d e+4 c f-3 c^2 f+6 d^2 e x+4 d f x+3 d^2 f x^2-8 f \log \left (\cos \left (\frac {1}{2} (c+d x)\right )+\sin \left (\frac {1}{2} (c+d x)\right )\right )-4 f \sin (c+d x)-d (e+f x) \sin (2 (c+d x))\right )\right )\right )}{8 a d^2 (1+\sin (c+d x))} \]

[In]

Integrate[((e + f*x)*Sin[c + d*x]^3)/(a + a*Sin[c + d*x]),x]

[Out]

((Cos[(c + d*x)/2] + Sin[(c + d*x)/2])*(Sin[(c + d*x)/2]*(8*d*(e + f*x)*Cos[c + d*x] - f*Cos[2*(c + d*x)] + 2*
(-8*d*e + 6*c*d*e + 4*c*f - 3*c^2*f + 6*d^2*e*x - 4*d*f*x + 3*d^2*f*x^2 - 8*f*Log[Cos[(c + d*x)/2] + Sin[(c +
d*x)/2]] - 4*f*Sin[c + d*x] - d*(e + f*x)*Sin[2*(c + d*x)])) + Cos[(c + d*x)/2]*(8*d*(e + f*x)*Cos[c + d*x] -
f*Cos[2*(c + d*x)] + 2*(6*c*d*e + 4*c*f - 3*c^2*f + 6*d^2*e*x + 4*d*f*x + 3*d^2*f*x^2 - 8*f*Log[Cos[(c + d*x)/
2] + Sin[(c + d*x)/2]] - 4*f*Sin[c + d*x] - d*(e + f*x)*Sin[2*(c + d*x)]))))/(8*a*d^2*(1 + Sin[c + d*x]))

Maple [C] (verified)

Result contains complex when optimal does not.

Time = 0.59 (sec) , antiderivative size = 187, normalized size of antiderivative = 1.18

method result size
risch \(\frac {3 f \,x^{2}}{4 a}+\frac {3 e x}{2 a}+\frac {\left (d x f +d e +i f \right ) {\mathrm e}^{i \left (d x +c \right )}}{2 a \,d^{2}}+\frac {\left (d x f +d e -i f \right ) {\mathrm e}^{-i \left (d x +c \right )}}{2 a \,d^{2}}+\frac {2 i f x}{a d}+\frac {2 i f c}{a \,d^{2}}+\frac {2 f x +2 e}{d a \left ({\mathrm e}^{i \left (d x +c \right )}+i\right )}-\frac {2 f \ln \left ({\mathrm e}^{i \left (d x +c \right )}+i\right )}{a \,d^{2}}-\frac {f \cos \left (2 d x +2 c \right )}{8 a \,d^{2}}-\frac {\left (f x +e \right ) \sin \left (2 d x +2 c \right )}{4 d a}\) \(187\)
parallelrisch \(\frac {16 f \left (\sin \left (\frac {d x}{2}+\frac {c}{2}\right )+\cos \left (\frac {d x}{2}+\frac {c}{2}\right )\right ) \ln \left (\sec ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )-32 \left (\sin \left (\frac {d x}{2}+\frac {c}{2}\right )+\cos \left (\frac {d x}{2}+\frac {c}{2}\right )\right ) f \ln \left (\tan \left (\frac {d x}{2}+\frac {c}{2}\right )+1\right )+\left (\left (12 d^{2} x^{2}+24 d x +26\right ) f +24 d^{2} e x +8 d e \right ) \cos \left (\frac {d x}{2}+\frac {c}{2}\right )+\left (\left (12 d^{2} x^{2}-24 d x +26\right ) f +24 d^{2} e x -40 d e \right ) \sin \left (\frac {d x}{2}+\frac {c}{2}\right )+\left (6 d x f +6 d e +7 f \right ) \cos \left (\frac {3 d x}{2}+\frac {3 c}{2}\right )+\left (2 d x f +2 d e -f \right ) \cos \left (\frac {5 d x}{2}+\frac {5 c}{2}\right )+\left (6 d x f +6 d e -7 f \right ) \sin \left (\frac {3 d x}{2}+\frac {3 c}{2}\right )-2 \sin \left (\frac {5 d x}{2}+\frac {5 c}{2}\right ) \left (d x f +d e +\frac {1}{2} f \right )}{16 a \,d^{2} \left (\sin \left (\frac {d x}{2}+\frac {c}{2}\right )+\cos \left (\frac {d x}{2}+\frac {c}{2}\right )\right )}\) \(266\)
default \(-\frac {-\frac {4 e \left (\frac {2}{\tan \left (\frac {d x}{2}+\frac {c}{2}\right )+1}+2 \arctan \left (\tan \left (\frac {d x}{2}+\frac {c}{2}\right )\right )\right )}{d}-\frac {2 f \,x^{2}+2 f \,x^{2} \tan \left (\frac {d x}{2}+\frac {c}{2}\right )+2 f \,x^{2} \left (\tan ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )+2 f \,x^{2} \left (\tan ^{3}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )+\frac {4 f x}{d}-\frac {4 f x \tan \left (\frac {d x}{2}+\frac {c}{2}\right )}{d}+\frac {4 f x \left (\tan ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{d}-\frac {4 f x \left (\tan ^{3}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{d}}{\left (\tan \left (\frac {d x}{2}+\frac {c}{2}\right )+1\right ) \left (1+\tan ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}+\frac {8 f \ln \left (\tan \left (\frac {d x}{2}+\frac {c}{2}\right )+1\right )}{d^{2}}-\frac {4 f \ln \left (1+\tan ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{d^{2}}+\frac {4 e \left (\frac {-\left (\tan ^{3}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )-2 \left (\tan ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )+\tan \left (\frac {d x}{2}+\frac {c}{2}\right )-2}{\left (1+\tan ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )^{2}}-\arctan \left (\tan \left (\frac {d x}{2}+\frac {c}{2}\right )\right )\right )}{d}+\frac {4 f \left (\left (d x +c \right ) \left (\frac {\cos \left (d x +c \right ) \sin \left (d x +c \right )}{2}+\frac {d x}{2}+\frac {c}{2}\right )-\frac {3 \left (d x +c \right )^{2}}{4}-\frac {\left (\sin ^{2}\left (d x +c \right )\right )}{4}+\sin \left (d x +c \right )-\cos \left (d x +c \right ) \left (d x +c \right )-c \left (\frac {\cos \left (d x +c \right ) \sin \left (d x +c \right )}{2}+\frac {d x}{2}+\frac {c}{2}\right )+\cos \left (d x +c \right ) c +\left (d x +c \right ) c \right )}{d^{2}}}{4 a}\) \(412\)
norman \(\frac {\frac {d e +2 f}{a \,d^{2}}-\frac {2 e \tan \left (\frac {d x}{2}+\frac {c}{2}\right )}{d a}+\frac {5 f \left (\tan ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{a \,d^{2}}+\frac {\left (-3 d e +2 f \right ) \left (\tan ^{7}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{a \,d^{2}}+\frac {\left (-6 d e +5 f \right ) \left (\tan ^{5}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{a \,d^{2}}+\frac {\left (-5 d e +3 f \right ) \left (\tan ^{3}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{a \,d^{2}}+\frac {\left (-d e +3 f \right ) \left (\tan ^{4}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{a \,d^{2}}+\frac {3 f \,x^{2}}{4 a}+\frac {3 f \,x^{2} \tan \left (\frac {d x}{2}+\frac {c}{2}\right )}{4 a}+\frac {9 f \,x^{2} \left (\tan ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{4 a}+\frac {9 f \,x^{2} \left (\tan ^{3}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{4 a}+\frac {9 f \,x^{2} \left (\tan ^{4}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{4 a}+\frac {9 f \,x^{2} \left (\tan ^{5}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{4 a}+\frac {3 f \,x^{2} \left (\tan ^{6}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{4 a}+\frac {3 f \,x^{2} \left (\tan ^{7}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{4 a}+\frac {\left (3 d e +4 f \right ) x}{2 d a}+\frac {\left (3 d e -4 f \right ) x \left (\tan ^{7}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{2 d a}+\frac {\left (3 d e -2 f \right ) x \tan \left (\frac {d x}{2}+\frac {c}{2}\right )}{2 d a}+\frac {3 \left (3 d e -2 f \right ) x \left (\tan ^{5}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{2 d a}+\frac {3 \left (3 d e +2 f \right ) x \left (\tan ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{2 d a}+\frac {\left (3 d e +2 f \right ) x \left (\tan ^{6}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{2 d a}+\frac {\left (9 d e -4 f \right ) x \left (\tan ^{3}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{2 d a}+\frac {\left (9 d e +4 f \right ) x \left (\tan ^{4}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{2 d a}}{\left (1+\tan ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )^{3} \left (\tan \left (\frac {d x}{2}+\frac {c}{2}\right )+1\right )}+\frac {f \ln \left (1+\tan ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{a \,d^{2}}-\frac {2 f \ln \left (\tan \left (\frac {d x}{2}+\frac {c}{2}\right )+1\right )}{a \,d^{2}}\) \(589\)

[In]

int((f*x+e)*sin(d*x+c)^3/(a+a*sin(d*x+c)),x,method=_RETURNVERBOSE)

[Out]

3/4*f*x^2/a+3/2*e*x/a+1/2*(d*x*f+I*f+d*e)/a/d^2*exp(I*(d*x+c))+1/2*(d*x*f-I*f+d*e)/a/d^2*exp(-I*(d*x+c))+2*I*f
/a/d*x+2*I*f/a/d^2*c+2*(f*x+e)/d/a/(exp(I*(d*x+c))+I)-2*f/a/d^2*ln(exp(I*(d*x+c))+I)-1/8*f/a/d^2*cos(2*d*x+2*c
)-1/4*(f*x+e)/d/a*sin(2*d*x+2*c)

Fricas [A] (verification not implemented)

none

Time = 0.27 (sec) , antiderivative size = 250, normalized size of antiderivative = 1.58 \[ \int \frac {(e+f x) \sin ^3(c+d x)}{a+a \sin (c+d x)} \, dx=\frac {6 \, d^{2} f x^{2} + 2 \, {\left (2 \, d f x + 2 \, d e - f\right )} \cos \left (d x + c\right )^{3} + 2 \, {\left (4 \, d f x + 4 \, d e + 3 \, f\right )} \cos \left (d x + c\right )^{2} + 8 \, d e + 4 \, {\left (3 \, d^{2} e + 2 \, d f\right )} x + {\left (6 \, d^{2} f x^{2} + 12 \, d e + 12 \, {\left (d^{2} e + d f\right )} x + f\right )} \cos \left (d x + c\right ) - 8 \, {\left (f \cos \left (d x + c\right ) + f \sin \left (d x + c\right ) + f\right )} \log \left (\sin \left (d x + c\right ) + 1\right ) + {\left (6 \, d^{2} f x^{2} - 2 \, {\left (2 \, d f x + 2 \, d e + f\right )} \cos \left (d x + c\right )^{2} - 8 \, d e + 4 \, {\left (3 \, d^{2} e - 2 \, d f\right )} x + 4 \, {\left (d f x + d e - 2 \, f\right )} \cos \left (d x + c\right ) - 7 \, f\right )} \sin \left (d x + c\right ) - 7 \, f}{8 \, {\left (a d^{2} \cos \left (d x + c\right ) + a d^{2} \sin \left (d x + c\right ) + a d^{2}\right )}} \]

[In]

integrate((f*x+e)*sin(d*x+c)^3/(a+a*sin(d*x+c)),x, algorithm="fricas")

[Out]

1/8*(6*d^2*f*x^2 + 2*(2*d*f*x + 2*d*e - f)*cos(d*x + c)^3 + 2*(4*d*f*x + 4*d*e + 3*f)*cos(d*x + c)^2 + 8*d*e +
 4*(3*d^2*e + 2*d*f)*x + (6*d^2*f*x^2 + 12*d*e + 12*(d^2*e + d*f)*x + f)*cos(d*x + c) - 8*(f*cos(d*x + c) + f*
sin(d*x + c) + f)*log(sin(d*x + c) + 1) + (6*d^2*f*x^2 - 2*(2*d*f*x + 2*d*e + f)*cos(d*x + c)^2 - 8*d*e + 4*(3
*d^2*e - 2*d*f)*x + 4*(d*f*x + d*e - 2*f)*cos(d*x + c) - 7*f)*sin(d*x + c) - 7*f)/(a*d^2*cos(d*x + c) + a*d^2*
sin(d*x + c) + a*d^2)

Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 4653 vs. \(2 (134) = 268\).

Time = 2.15 (sec) , antiderivative size = 4653, normalized size of antiderivative = 29.45 \[ \int \frac {(e+f x) \sin ^3(c+d x)}{a+a \sin (c+d x)} \, dx=\text {Too large to display} \]

[In]

integrate((f*x+e)*sin(d*x+c)**3/(a+a*sin(d*x+c)),x)

[Out]

Piecewise((6*d**2*e*x*tan(c/2 + d*x/2)**5/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d
**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 6*d**2*e*x*ta
n(c/2 + d*x/2)**4/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3
+ 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 12*d**2*e*x*tan(c/2 + d*x/2)**3/(4*a*
d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*
x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 12*d**2*e*x*tan(c/2 + d*x/2)**2/(4*a*d**2*tan(c/2 + d*x/2)**
5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(
c/2 + d*x/2) + 4*a*d**2) + 6*d**2*e*x*tan(c/2 + d*x/2)/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/
2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) +
6*d**2*e*x/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d
**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 3*d**2*f*x**2*tan(c/2 + d*x/2)**5/(4*a*d**2*
tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)*
*2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 3*d**2*f*x**2*tan(c/2 + d*x/2)**4/(4*a*d**2*tan(c/2 + d*x/2)**5 +
 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2
 + d*x/2) + 4*a*d**2) + 6*d**2*f*x**2*tan(c/2 + d*x/2)**3/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d
*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2)
 + 6*d**2*f*x**2*tan(c/2 + d*x/2)**2/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*t
an(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 3*d**2*f*x**2*tan(
c/2 + d*x/2)/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a
*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 3*d**2*f*x**2/(4*a*d**2*tan(c/2 + d*x/2)**
5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(
c/2 + d*x/2) + 4*a*d**2) + 12*d*e*tan(c/2 + d*x/2)**4/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2
)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 1
2*d*e*tan(c/2 + d*x/2)**3/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*
x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 20*d*e*tan(c/2 + d*x/2)**2/(4
*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 +
 d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 4*d*e*tan(c/2 + d*x/2)/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*
a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 +
d*x/2) + 4*a*d**2) + 16*d*e/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 +
d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) - 8*d*f*x*tan(c/2 + d*x/2)**5
/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/
2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 4*d*f*x*tan(c/2 + d*x/2)**4/(4*a*d**2*tan(c/2 + d*x/2)
**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*ta
n(c/2 + d*x/2) + 4*a*d**2) - 4*d*f*x*tan(c/2 + d*x/2)**3/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*
x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2)
+ 4*d*f*x*tan(c/2 + d*x/2)**2/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2
+ d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) - 4*d*f*x*tan(c/2 + d*x/2)/
(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2
 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 8*d*f*x/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/
2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*
d**2) - 8*f*log(tan(c/2 + d*x/2) + 1)*tan(c/2 + d*x/2)**5/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d
*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2)
 - 8*f*log(tan(c/2 + d*x/2) + 1)*tan(c/2 + d*x/2)**4/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)
**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) - 16
*f*log(tan(c/2 + d*x/2) + 1)*tan(c/2 + d*x/2)**3/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4
+ 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) - 16*f*l
og(tan(c/2 + d*x/2) + 1)*tan(c/2 + d*x/2)**2/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*
a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) - 8*f*log(ta
n(c/2 + d*x/2) + 1)*tan(c/2 + d*x/2)/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*t
an(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) - 8*f*log(tan(c/2 +
d*x/2) + 1)/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*
d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 4*f*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 +
d*x/2)**5/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d*
*2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 4*f*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*
x/2)**4/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2
*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 8*f*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/
2)**3/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*t
an(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 8*f*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)
**2/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan
(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 4*f*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)/(
4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2
+ d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) + 4*f*log(tan(c/2 + d*x/2)**2 + 1)/(4*a*d**2*tan(c/2 + d*x
/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2
*tan(c/2 + d*x/2) + 4*a*d**2) - 8*f*tan(c/2 + d*x/2)**4/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x
/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) -
 4*f*tan(c/2 + d*x/2)**3/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x
/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) - 4*f*tan(c/2 + d*x/2)**2/(4*a*d
**2*tan(c/2 + d*x/2)**5 + 4*a*d**2*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x
/2)**2 + 4*a*d**2*tan(c/2 + d*x/2) + 4*a*d**2) - 8*f*tan(c/2 + d*x/2)/(4*a*d**2*tan(c/2 + d*x/2)**5 + 4*a*d**2
*tan(c/2 + d*x/2)**4 + 8*a*d**2*tan(c/2 + d*x/2)**3 + 8*a*d**2*tan(c/2 + d*x/2)**2 + 4*a*d**2*tan(c/2 + d*x/2)
 + 4*a*d**2), Ne(d, 0)), ((e*x + f*x**2/2)*sin(c)**3/(a*sin(c) + a), True))

Maxima [F(-2)]

Exception generated. \[ \int \frac {(e+f x) \sin ^3(c+d x)}{a+a \sin (c+d x)} \, dx=\text {Exception raised: RuntimeError} \]

[In]

integrate((f*x+e)*sin(d*x+c)^3/(a+a*sin(d*x+c)),x, algorithm="maxima")

[Out]

Exception raised: RuntimeError >> ECL says: THROW: The catch RAT-ERR is undefined.

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 7564 vs. \(2 (138) = 276\).

Time = 0.98 (sec) , antiderivative size = 7564, normalized size of antiderivative = 47.87 \[ \int \frac {(e+f x) \sin ^3(c+d x)}{a+a \sin (c+d x)} \, dx=\text {Too large to display} \]

[In]

integrate((f*x+e)*sin(d*x+c)^3/(a+a*sin(d*x+c)),x, algorithm="giac")

[Out]

1/8*(6*d^2*f*x^2*tan(1/2*d*x)^5*tan(1/2*c)^5 - 6*d^2*f*x^2*tan(1/2*d*x)^5*tan(1/2*c)^4 - 6*d^2*f*x^2*tan(1/2*d
*x)^4*tan(1/2*c)^5 + 12*d^2*e*x*tan(1/2*d*x)^5*tan(1/2*c)^5 + 12*d^2*f*x^2*tan(1/2*d*x)^5*tan(1/2*c)^3 - 6*d^2
*f*x^2*tan(1/2*d*x)^4*tan(1/2*c)^4 - 12*d^2*e*x*tan(1/2*d*x)^5*tan(1/2*c)^4 + 12*d^2*f*x^2*tan(1/2*d*x)^3*tan(
1/2*c)^5 - 12*d^2*e*x*tan(1/2*d*x)^4*tan(1/2*c)^5 + 16*d*f*x*tan(1/2*d*x)^5*tan(1/2*c)^5 - 12*d^2*f*x^2*tan(1/
2*d*x)^5*tan(1/2*c)^2 - 12*d^2*f*x^2*tan(1/2*d*x)^4*tan(1/2*c)^3 + 24*d^2*e*x*tan(1/2*d*x)^5*tan(1/2*c)^3 - 12
*d^2*f*x^2*tan(1/2*d*x)^3*tan(1/2*c)^4 - 12*d^2*e*x*tan(1/2*d*x)^4*tan(1/2*c)^4 + 8*d*f*x*tan(1/2*d*x)^5*tan(1
/2*c)^4 - 12*d^2*f*x^2*tan(1/2*d*x)^2*tan(1/2*c)^5 + 24*d^2*e*x*tan(1/2*d*x)^3*tan(1/2*c)^5 + 8*d*f*x*tan(1/2*
d*x)^4*tan(1/2*c)^5 + 16*d*e*tan(1/2*d*x)^5*tan(1/2*c)^5 - 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*
d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c
) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c)^5 + 6*d^2*
f*x^2*tan(1/2*d*x)^5*tan(1/2*c) - 12*d^2*f*x^2*tan(1/2*d*x)^4*tan(1/2*c)^2 - 24*d^2*e*x*tan(1/2*d*x)^5*tan(1/2
*c)^2 + 24*d^2*f*x^2*tan(1/2*d*x)^3*tan(1/2*c)^3 - 24*d^2*e*x*tan(1/2*d*x)^4*tan(1/2*c)^3 + 8*d*f*x*tan(1/2*d*
x)^5*tan(1/2*c)^3 - 12*d^2*f*x^2*tan(1/2*d*x)^2*tan(1/2*c)^4 - 24*d^2*e*x*tan(1/2*d*x)^3*tan(1/2*c)^4 - 64*d*f
*x*tan(1/2*d*x)^4*tan(1/2*c)^4 + 8*d*e*tan(1/2*d*x)^5*tan(1/2*c)^4 + 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 -
2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2
*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c)^
4 + 6*d^2*f*x^2*tan(1/2*d*x)*tan(1/2*c)^5 - 24*d^2*e*x*tan(1/2*d*x)^2*tan(1/2*c)^5 + 8*d*f*x*tan(1/2*d*x)^3*ta
n(1/2*c)^5 + 8*d*e*tan(1/2*d*x)^4*tan(1/2*c)^5 + 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan
(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(ta
n(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4*tan(1/2*c)^5 - f*tan(1/2*d*x)^5
*tan(1/2*c)^5 - 6*d^2*f*x^2*tan(1/2*d*x)^5 - 6*d^2*f*x^2*tan(1/2*d*x)^4*tan(1/2*c) + 12*d^2*e*x*tan(1/2*d*x)^5
*tan(1/2*c) - 24*d^2*f*x^2*tan(1/2*d*x)^3*tan(1/2*c)^2 - 24*d^2*e*x*tan(1/2*d*x)^4*tan(1/2*c)^2 + 8*d*f*x*tan(
1/2*d*x)^5*tan(1/2*c)^2 - 24*d^2*f*x^2*tan(1/2*d*x)^2*tan(1/2*c)^3 + 48*d^2*e*x*tan(1/2*d*x)^3*tan(1/2*c)^3 -
8*d*f*x*tan(1/2*d*x)^4*tan(1/2*c)^3 + 8*d*e*tan(1/2*d*x)^5*tan(1/2*c)^3 - 16*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c
)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*
x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1
/2*c)^3 - 6*d^2*f*x^2*tan(1/2*d*x)*tan(1/2*c)^4 - 24*d^2*e*x*tan(1/2*d*x)^2*tan(1/2*c)^4 - 8*d*f*x*tan(1/2*d*x
)^3*tan(1/2*c)^4 - 64*d*e*tan(1/2*d*x)^4*tan(1/2*c)^4 + 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x
)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) +
 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4*tan(1/2*c)^4 + 17*f*tan(
1/2*d*x)^5*tan(1/2*c)^4 - 6*d^2*f*x^2*tan(1/2*c)^5 + 12*d^2*e*x*tan(1/2*d*x)*tan(1/2*c)^5 + 8*d*f*x*tan(1/2*d*
x)^2*tan(1/2*c)^5 + 8*d*e*tan(1/2*d*x)^3*tan(1/2*c)^5 - 16*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*
x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c)
+ 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c)^5 + 17*f*tan
(1/2*d*x)^4*tan(1/2*c)^5 - 6*d^2*f*x^2*tan(1/2*d*x)^4 - 12*d^2*e*x*tan(1/2*d*x)^5 + 12*d^2*f*x^2*tan(1/2*d*x)^
3*tan(1/2*c) - 12*d^2*e*x*tan(1/2*d*x)^4*tan(1/2*c) + 8*d*f*x*tan(1/2*d*x)^5*tan(1/2*c) - 24*d^2*f*x^2*tan(1/2
*d*x)^2*tan(1/2*c)^2 - 48*d^2*e*x*tan(1/2*d*x)^3*tan(1/2*c)^2 + 8*d*f*x*tan(1/2*d*x)^4*tan(1/2*c)^2 + 8*d*e*ta
n(1/2*d*x)^5*tan(1/2*c)^2 + 16*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*
d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2
*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c)^2 + 12*d^2*f*x^2*tan(1/2*d*x)*tan(1/2*c)
^3 - 48*d^2*e*x*tan(1/2*d*x)^2*tan(1/2*c)^3 + 160*d*f*x*tan(1/2*d*x)^3*tan(1/2*c)^3 - 8*d*e*tan(1/2*d*x)^4*tan
(1/2*c)^3 + 16*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^
2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*
d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4*tan(1/2*c)^3 - 10*f*tan(1/2*d*x)^5*tan(1/2*c)^3 - 6*d^2*f*x^2*tan(1
/2*c)^4 - 12*d^2*e*x*tan(1/2*d*x)*tan(1/2*c)^4 + 8*d*f*x*tan(1/2*d*x)^2*tan(1/2*c)^4 - 8*d*e*tan(1/2*d*x)^3*ta
n(1/2*c)^4 + 16*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)
^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2
*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c)^4 - 15*f*tan(1/2*d*x)^4*tan(1/2*c)^4 - 12*d^2*e*x*tan(1
/2*c)^5 + 8*d*f*x*tan(1/2*d*x)*tan(1/2*c)^5 + 8*d*e*tan(1/2*d*x)^2*tan(1/2*c)^5 + 16*f*log(2*(tan(1/2*d*x)^2*t
an(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*ta
n(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)
^2*tan(1/2*c)^5 - 10*f*tan(1/2*d*x)^3*tan(1/2*c)^5 - 12*d^2*f*x^2*tan(1/2*d*x)^3 - 12*d^2*e*x*tan(1/2*d*x)^4 +
 16*d*f*x*tan(1/2*d*x)^5 - 12*d^2*f*x^2*tan(1/2*d*x)^2*tan(1/2*c) + 24*d^2*e*x*tan(1/2*d*x)^3*tan(1/2*c) + 64*
d*f*x*tan(1/2*d*x)^4*tan(1/2*c) + 8*d*e*tan(1/2*d*x)^5*tan(1/2*c) - 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2
*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*
tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c) -
 12*d^2*f*x^2*tan(1/2*d*x)*tan(1/2*c)^2 - 48*d^2*e*x*tan(1/2*d*x)^2*tan(1/2*c)^2 + 160*d*f*x*tan(1/2*d*x)^3*ta
n(1/2*c)^2 + 8*d*e*tan(1/2*d*x)^4*tan(1/2*c)^2 + 16*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*ta
n(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(t
an(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4*tan(1/2*c)^2 + 10*f*tan(1/2*d*
x)^5*tan(1/2*c)^2 - 12*d^2*f*x^2*tan(1/2*c)^3 + 24*d^2*e*x*tan(1/2*d*x)*tan(1/2*c)^3 + 160*d*f*x*tan(1/2*d*x)^
2*tan(1/2*c)^3 + 160*d*e*tan(1/2*d*x)^3*tan(1/2*c)^3 - 32*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x
)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) +
 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c)^3 - 38*f*tan(
1/2*d*x)^4*tan(1/2*c)^3 - 12*d^2*e*x*tan(1/2*c)^4 + 64*d*f*x*tan(1/2*d*x)*tan(1/2*c)^4 + 8*d*e*tan(1/2*d*x)^2*
tan(1/2*c)^4 + 16*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*
c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1
/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^4 - 38*f*tan(1/2*d*x)^3*tan(1/2*c)^4 + 16*d*f*x*tan(1
/2*c)^5 + 8*d*e*tan(1/2*d*x)*tan(1/2*c)^5 - 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*
c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2
*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*tan(1/2*c)^5 + 10*f*tan(1/2*d*x)^2*tan
(1/2*c)^5 - 12*d^2*f*x^2*tan(1/2*d*x)^2 - 24*d^2*e*x*tan(1/2*d*x)^3 - 8*d*f*x*tan(1/2*d*x)^4 + 16*d*e*tan(1/2*
d*x)^5 + 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 +
tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)
^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5 + 6*d^2*f*x^2*tan(1/2*d*x)*tan(1/2*c) - 24*d^2*e*x*tan(1/2*d*x)^2*tan(1
/2*c) - 8*d*f*x*tan(1/2*d*x)^3*tan(1/2*c) + 64*d*e*tan(1/2*d*x)^4*tan(1/2*c) + 8*f*log(2*(tan(1/2*d*x)^2*tan(1
/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/
2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4*t
an(1/2*c) - 17*f*tan(1/2*d*x)^5*tan(1/2*c) - 12*d^2*f*x^2*tan(1/2*c)^2 - 24*d^2*e*x*tan(1/2*d*x)*tan(1/2*c)^2
- 160*d*f*x*tan(1/2*d*x)^2*tan(1/2*c)^2 + 160*d*e*tan(1/2*d*x)^3*tan(1/2*c)^2 + 32*f*log(2*(tan(1/2*d*x)^2*tan
(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(
1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3
*tan(1/2*c)^2 - 38*f*tan(1/2*d*x)^4*tan(1/2*c)^2 - 24*d^2*e*x*tan(1/2*c)^3 - 8*d*f*x*tan(1/2*d*x)*tan(1/2*c)^3
 + 160*d*e*tan(1/2*d*x)^2*tan(1/2*c)^3 + 32*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c)
 - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d
*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^3 - 52*f*tan(1/2*d*x)^3*tan
(1/2*c)^3 - 8*d*f*x*tan(1/2*c)^4 + 64*d*e*tan(1/2*d*x)*tan(1/2*c)^4 + 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 -
 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) +
2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*tan(1/2*c)^4
 - 38*f*tan(1/2*d*x)^2*tan(1/2*c)^4 + 16*d*e*tan(1/2*c)^5 + 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2
*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*
c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^5 - 17*f*tan(1/2*d*x)*ta
n(1/2*c)^5 - 6*d^2*f*x^2*tan(1/2*d*x) - 24*d^2*e*x*tan(1/2*d*x)^2 + 8*d*f*x*tan(1/2*d*x)^3 - 8*d*e*tan(1/2*d*x
)^4 + 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan
(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2
+ tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4 + f*tan(1/2*d*x)^5 - 6*d^2*f*x^2*tan(1/2*c) + 12*d^2*e*x*tan(1/2*d*x)*tan(
1/2*c) - 8*d*f*x*tan(1/2*d*x)^2*tan(1/2*c) - 8*d*e*tan(1/2*d*x)^3*tan(1/2*c) - 16*f*log(2*(tan(1/2*d*x)^2*tan(
1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1
/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*
tan(1/2*c) - 15*f*tan(1/2*d*x)^4*tan(1/2*c) - 24*d^2*e*x*tan(1/2*c)^2 - 8*d*f*x*tan(1/2*d*x)*tan(1/2*c)^2 - 16
0*d*e*tan(1/2*d*x)^2*tan(1/2*c)^2 + 32*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*
tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2
*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^2 + 52*f*tan(1/2*d*x)^3*tan(1/2*
c)^2 + 8*d*f*x*tan(1/2*c)^3 - 8*d*e*tan(1/2*d*x)*tan(1/2*c)^3 - 16*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*ta
n(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan
(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*tan(1/2*c)^3 + 52
*f*tan(1/2*d*x)^2*tan(1/2*c)^3 - 8*d*e*tan(1/2*c)^4 + 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^
2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1
)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^4 - 15*f*tan(1/2*d*x)*tan(1/2*
c)^4 + f*tan(1/2*c)^5 - 6*d^2*f*x^2 - 12*d^2*e*x*tan(1/2*d*x) - 8*d*f*x*tan(1/2*d*x)^2 + 8*d*e*tan(1/2*d*x)^3
+ 16*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/
2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + t
an(1/2*c)^2 + 1))*tan(1/2*d*x)^3 + 17*f*tan(1/2*d*x)^4 - 12*d^2*e*x*tan(1/2*c) + 64*d*f*x*tan(1/2*d*x)*tan(1/2
*c) - 8*d*e*tan(1/2*d*x)^2*tan(1/2*c) + 16*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c)
- 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*
x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c) + 38*f*tan(1/2*d*x)^3*tan(1/
2*c) - 8*d*f*x*tan(1/2*c)^2 - 8*d*e*tan(1/2*d*x)*tan(1/2*c)^2 + 16*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*ta
n(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan
(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*tan(1/2*c)^2 + 52
*f*tan(1/2*d*x)^2*tan(1/2*c)^2 + 8*d*e*tan(1/2*c)^3 + 16*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)
^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) +
1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^3 + 38*f*tan(1/2*d*x)*tan(1/2
*c)^3 + 17*f*tan(1/2*c)^4 - 12*d^2*e*x + 8*d*f*x*tan(1/2*d*x) - 8*d*e*tan(1/2*d*x)^2 + 16*f*log(2*(tan(1/2*d*x
)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 +
 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2
*d*x)^2 + 10*f*tan(1/2*d*x)^3 + 8*d*f*x*tan(1/2*c) + 64*d*e*tan(1/2*d*x)*tan(1/2*c) - 8*f*log(2*(tan(1/2*d*x)^
2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2
*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d
*x)*tan(1/2*c) - 38*f*tan(1/2*d*x)^2*tan(1/2*c) - 8*d*e*tan(1/2*c)^2 + 16*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2
 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x)
+ 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^2 - 38*f*tan
(1/2*d*x)*tan(1/2*c)^2 + 10*f*tan(1/2*c)^3 - 16*d*f*x + 8*d*e*tan(1/2*d*x) + 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2
*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*
d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x) + 10*
f*tan(1/2*d*x)^2 + 8*d*e*tan(1/2*c) + 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2
*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^
2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c) + 15*f*tan(1/2*d*x)*tan(1/2*c) + 10*f*tan(1/2*
c)^2 - 16*d*e + 8*f*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*
c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1
/2*d*x)^2 + tan(1/2*c)^2 + 1)) + 17*f*tan(1/2*d*x) + 17*f*tan(1/2*c) + f)/(a*d^2*tan(1/2*d*x)^5*tan(1/2*c)^5 -
 a*d^2*tan(1/2*d*x)^5*tan(1/2*c)^4 - a*d^2*tan(1/2*d*x)^4*tan(1/2*c)^5 + 2*a*d^2*tan(1/2*d*x)^5*tan(1/2*c)^3 -
 a*d^2*tan(1/2*d*x)^4*tan(1/2*c)^4 + 2*a*d^2*tan(1/2*d*x)^3*tan(1/2*c)^5 - 2*a*d^2*tan(1/2*d*x)^5*tan(1/2*c)^2
 - 2*a*d^2*tan(1/2*d*x)^4*tan(1/2*c)^3 - 2*a*d^2*tan(1/2*d*x)^3*tan(1/2*c)^4 - 2*a*d^2*tan(1/2*d*x)^2*tan(1/2*
c)^5 + a*d^2*tan(1/2*d*x)^5*tan(1/2*c) - 2*a*d^2*tan(1/2*d*x)^4*tan(1/2*c)^2 + 4*a*d^2*tan(1/2*d*x)^3*tan(1/2*
c)^3 - 2*a*d^2*tan(1/2*d*x)^2*tan(1/2*c)^4 + a*d^2*tan(1/2*d*x)*tan(1/2*c)^5 - a*d^2*tan(1/2*d*x)^5 - a*d^2*ta
n(1/2*d*x)^4*tan(1/2*c) - 4*a*d^2*tan(1/2*d*x)^3*tan(1/2*c)^2 - 4*a*d^2*tan(1/2*d*x)^2*tan(1/2*c)^3 - a*d^2*ta
n(1/2*d*x)*tan(1/2*c)^4 - a*d^2*tan(1/2*c)^5 - a*d^2*tan(1/2*d*x)^4 + 2*a*d^2*tan(1/2*d*x)^3*tan(1/2*c) - 4*a*
d^2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*a*d^2*tan(1/2*d*x)*tan(1/2*c)^3 - a*d^2*tan(1/2*c)^4 - 2*a*d^2*tan(1/2*d*x
)^3 - 2*a*d^2*tan(1/2*d*x)^2*tan(1/2*c) - 2*a*d^2*tan(1/2*d*x)*tan(1/2*c)^2 - 2*a*d^2*tan(1/2*c)^3 - 2*a*d^2*t
an(1/2*d*x)^2 + a*d^2*tan(1/2*d*x)*tan(1/2*c) - 2*a*d^2*tan(1/2*c)^2 - a*d^2*tan(1/2*d*x) - a*d^2*tan(1/2*c) -
 a*d^2)

Mupad [B] (verification not implemented)

Time = 1.90 (sec) , antiderivative size = 246, normalized size of antiderivative = 1.56 \[ \int \frac {(e+f x) \sin ^3(c+d x)}{a+a \sin (c+d x)} \, dx={\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}\,\left (\frac {d\,e+f\,1{}\mathrm {i}}{2\,a\,d^2}+\frac {f\,x}{2\,a\,d}\right )-{\mathrm {e}}^{-c\,1{}\mathrm {i}-d\,x\,1{}\mathrm {i}}\,\left (\frac {-d\,e+f\,1{}\mathrm {i}}{2\,a\,d^2}-\frac {f\,x}{2\,a\,d}\right )+{\mathrm {e}}^{-c\,2{}\mathrm {i}-d\,x\,2{}\mathrm {i}}\,\left (\frac {\left (-2\,d\,e+f\,1{}\mathrm {i}\right )\,1{}\mathrm {i}}{16\,a\,d^2}-\frac {f\,x\,1{}\mathrm {i}}{8\,a\,d}\right )+{\mathrm {e}}^{c\,2{}\mathrm {i}+d\,x\,2{}\mathrm {i}}\,\left (\frac {\left (2\,d\,e+f\,1{}\mathrm {i}\right )\,1{}\mathrm {i}}{16\,a\,d^2}+\frac {f\,x\,1{}\mathrm {i}}{8\,a\,d}\right )+\frac {3\,f\,x^2}{4\,a}-\frac {2\,f\,\ln \left ({\mathrm {e}}^{c\,1{}\mathrm {i}}\,{\mathrm {e}}^{d\,x\,1{}\mathrm {i}}+1{}\mathrm {i}\right )}{a\,d^2}+\frac {2\,\left (e+f\,x\right )}{a\,d\,\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}+1{}\mathrm {i}\right )}+\frac {x\,\left (3\,d\,e+f\,4{}\mathrm {i}\right )}{2\,a\,d} \]

[In]

int((sin(c + d*x)^3*(e + f*x))/(a + a*sin(c + d*x)),x)

[Out]

exp(c*1i + d*x*1i)*((f*1i + d*e)/(2*a*d^2) + (f*x)/(2*a*d)) - exp(- c*1i - d*x*1i)*((f*1i - d*e)/(2*a*d^2) - (
f*x)/(2*a*d)) + exp(- c*2i - d*x*2i)*(((f*1i - 2*d*e)*1i)/(16*a*d^2) - (f*x*1i)/(8*a*d)) + exp(c*2i + d*x*2i)*
(((f*1i + 2*d*e)*1i)/(16*a*d^2) + (f*x*1i)/(8*a*d)) + (3*f*x^2)/(4*a) - (2*f*log(exp(c*1i)*exp(d*x*1i) + 1i))/
(a*d^2) + (2*(e + f*x))/(a*d*(exp(c*1i + d*x*1i) + 1i)) + (x*(f*4i + 3*d*e))/(2*a*d)